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Interatomic coulombic decay (ICD) is a relaxation process induced by electronic correlation. In 
this work we study the ICD process in a two coupled Quantum wells (QWs) nano-structure. We 
study a simple one-dimensional effective potential using experimental parameters of the semicon¬ 
ductor QW layers i.e. using the single band effective-mass approximation . In our calculations we 
consider the discontinuity of the effective mass of the electron in each of the QW layers. We control 
the ICD lifetime by changing the distance between the two wells. The expected overall trend is 
a decrease of ICD lifetime with a decrease in the distance between the wells. We show that the 
distance can be tuned such that the emitted ICD electron is trapped in a meta-stable state in the 
continuum i.e. a one electron resonance state. This causes the life time of the ICD to be an order 
of magnitude smaller even in very long distances, and improves the efficiency of the ICD. For the 
ICD to be dominant decay mechanism it must prevail over all other possible competitive decay 
processes. We have found that the lifetime of the ICD is on the timescale of picoseconds. Therefore, 
based on our results we can design an experiment that will observe the ICD phenomenon in QWs 
nano-structure for the first time. This work can lead to designing a wavelength sensitive detector 
which is efficient even in low intensities. 


I. INTRODUCTION 

Interatomic/intermolecular coulombic decay (ICD) is 
a very efficient and fast electron relaxation process re¬ 
lying on the correlation between electrons. Such pro¬ 
cess occurs by passing the excess excitation energy of 
one electron to another electron in a neighbouring atom 
or molecule, resulting in the ionization of this electron. 
ICD was first proposed by Cederbaum and co-workers in 
hydrogen bonded molecular clustersflj. Past studies fo¬ 
cused on weakly bounded systems such as van der Waals 
clusters and weakly bounded dimers like the he¬ 

lium dimer which is the most weakly bound system in 
nature Q. 

ICD was observed experimentally on nobel gas clus¬ 
ters and dimers like Ne, Ar and He, and in hydrogen- 
bonded systems like water molecule dimers [6}-|12|. All 
past studies show that the energy transfer in ICD through 
electron correlation happens also at extremely long dis¬ 
tances. This can occur due to the fact that the ionized 
electron has a long De-Broglie wavelength so its wave- 
function couples to the bound states involved in the pro¬ 
cess. 

There are few ways to trigger the ICD process. It could 
be produced directly from photo-ionization of an inner- 
valence electron j|] or as a result of multistage process 
such as photo-ionization followed by Auger ionization 
i EH- Recently another multistage resonant-Auger- 
driven ICD was proposed which does not involve photo¬ 
ionization of inner-valence electron but requires just the 
excitation of this electron to an unoccupied orbital |l4j . 
This special process can yield a very high sensitivity to 
the location and energies of the ICD electrons, and can 
be studied in big molecules such as proteins and DNA. 
It was shown both theoretically and experimentally that 
the inter-atomic decay rate is strongly dependent dis¬ 


tance of neighboring atoms[2, 1 , 5 ). In this sense, the 
ICD is more efficient as the distance between the atoms 
is smaller. Moreover it was shown that the ICD lifetime 
decreases as the number of neighbors increases 0] (l6| . 

Recently it was shown that coupled quantum dots can 
undergo ICD. Quantum dots (QDs) are solid structures 
composed of semiconductors which confines electrons in 
three dimensions and as such they serve as an artificial 
atoms 0B The ICD process in QDs was proven to be 
very efficient, in comparison to other decay mechanisms 
exists in the QD, having lifetime of picosecond and less 
0-111. It was shown that the ICD lifetime in QDs 
grows with the distance between the dots. Quantum wells 
(QWs) and QDs are widely used in optoelectronic devices 
such as laser diodes and photo-detectors. Compared to 
QDs, QWs are easier to grow and control its dimensions. 

In this work, we study the ICD decay process in a nano¬ 
structure composed of two coupled QWs with different 
widths, the ICD works in the same manner like in molec¬ 
ular clusters and QDs. An excited electron in one well 
passes its extra energy to the electron in the neighbor¬ 
ing well which is ionized. Excited electrons in QWs have 
many relaxation pathways, such as spontaneous photon 
emission, and interaction with phonons. These processes 
are competing with the ICD in our system. For the ICD 
to be the dominant decay process it must have lifetime 
on the same timescale of the shortest decay process in the 
system. Inter-band spontaneous photon emission is not 
the most efficient decay process in QW’s, and has life¬ 
time on the time scale of nsec [01 ■ Scattering of electron 
with longitudinal optic (LO) phonons i.e. inter-subband 
relaxation through vibrational modes of the lattice is the 
most efficient competing relaxation process with lifetime 
in the timescale of picosecond [0-f26j. 

In our calculation we use an effective one-dimensional 
potential which takes into account only one conduction 
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band of the QWs nano-structure i.e. we use the single 
band effective mass approximation. We compared the 
bound state energies and wave functions of this effective 
one electron hamiltonian with the bound states calcu¬ 
lated from the k ■ p method [28|- 32]. The k ■ p method 


takes into account contribution of eight bands from the 
conduction and valence bands. Each QW produces a 
rectangular potential well in the conduction band, which 
its depth is calculated from experimental measurements 
of the band gaps of each layer in the nano-structure 


We also take into account the change in effective mass 
of the electron. This effective mass is not constant over 
the entire space but depends on the semiconductor lay¬ 
ers that produces the QWs nano-structure. Here, for the 
first time, the real band structure and the discontinu¬ 
ity of the effective mass are taken into account in the 
ICD calculation. Furthermore, our QW nano-structure 
can be grown in the lab in an easy manner compared 
to QD systems. In this context the QWs widths and 
the distance between the wells can be controlled during 
the growth process. Due to the fact that the parameters 
of the QW nano-structure in the calculation are taken 
from real semiconductor properties, one can carry out an 
experiment based on the results of this paper. A brief 
description of such experiment will be given below. 

We study the ICD process here as a two-electron reso¬ 
nance function which has a finite lifetime. We calculate 
the life-time of the ICD resonance states at different dis¬ 
tances between the neighboring wells. We report an un¬ 
expected result, showing that although the wells are far 
apart from each other the lifetime of the ICD resonance 
is an order of magnitude shorter compared with the ex¬ 
pected value. This unexpected short lifetime enhances 
the efficiency of the ICD process over other competitive 
relaxation processes even when the wells are far apart and 
tunneling is unlikely to occur in these distances. We ex¬ 
plain the result using the one electron resonance energies 
of the ionized electron (38j| . This one-electron resonance 
originates from the rectangular shape of the two QWs, 
and introduces a large density of continuum states which 
makes the two-electron ICD process efficient and fast. 

The paper is organized as follows: In SecHT] we intro¬ 
duce a schematic representation of the ICD in the dou¬ 
ble QW nano-structure. We then move to the methods 
of calculating the ICD lifetime and a proper representa¬ 
tion of the hamiltonian in Sec lIIIl In Sec lIVI we present 
our finding and explain our results before concluding in 
Sec|V] 


II. SCHEMATIC REPRESENTATION OF THE 
ICD PROCESS IN THE DOUBLE QUANTUM 
WELL NANO-STRUCTURE 


The ICD process we study in this work is based on 
two electrons in two coupled QWs. A schematic presen¬ 
tation of one QW nano-structure is shown in FigJT] this 


QW consists of semiconductors and can be grown in the 
lab. By attaching different layers of semiconductor mate¬ 
rials such that a semiconductor with a smaller band gap 
is sandwiched between two semiconductor with a larger 
band gap a QW structure is created (see Fig[l]). Due to 
the difference in the band gaps of each material in the dif¬ 
ferent layers a well structure is formed in the conduction 
band, and a barrier is formed in the valence band. 

In our system the semiconductor with the smaller band 
gap is Ino. 53 Gao. 47 As. This semiconductor has band gap 
of Eg = 0.74eV and the effective mass of the electron 
in this layer is = 0.045m e . The material with the 
larger band gap is Ino. 52 Alo. 48 As, which has a band gap of 
Eg = 1.45eV and the effective mass of the electron in this 
layer is m h e j = 0.075m e [Hj. Both materials are lattice 
matched to InP, this eliminates all the strain effects in 
the system. 
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FIG. 1 . QW nano-structure composed of two different semi¬ 
conductor materials thin layers, we have chosen the growth 
direction to be in along the x axis. The semiconductor with 
the smaller band gap is sandwiched between two semicon¬ 
ductors with a larger band gap. This spacial structure trans¬ 
formed a well in the conduction band, and a barrier in the va¬ 
lence band of the nano-structure. The semiconductor with the 
smaller band gap is In2,2=o.53Ga y 2=o.47As and has a band gap 
of Eg = 0 . 74 eV. The semiconductor with the larger band gap 
is Inxi=o.52Alyi—o.48As and has a band gap of E g = 1 . 45 eV. 
The effective mass of the electron in Ino.53Gao.47As and 
Ino 52AI0 48As is m™ f = 0 . 045 m e , m b ef = 0.075 m e respectively 
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The nano-structure we study of two coupled QWs is 
shown schematically in FigJ5] The left QW is wider than 
the right QW and supports two bound states , while the 
right QW supports only one bound state. The ICD pro¬ 
cess in this system is based on the correlation between 
two electrons, each one in a different QW. To initiate 
the ICD process we excite two electrons from the valence 
band to the conduction band. We propose an experi¬ 
ment to measure the ICD process using two pulsed lasers 
at low temperature. In Fig Q](a)| we see a schematic repre¬ 
sentation of the conduction band and valence band of the 
coupled QWs nano-structure. The first laser with a fre¬ 
quency of uj\ matches the interband transition between 
the first hole bound state and the first excited electronic 
bound state located in the wider well (see Fig QKa)[ ). The 
second laser with a frequency of ui 2 matches the inter¬ 
band transition between the first hole bound state and 
the electronic bound state located in the narrow well (see 
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Fig Qta)[ ) . By aligning the two lasers to excite the same 
spot (perpendicular to the growth direction) on the sam¬ 
ple, we can treat the problem as one dimensional. This 
first step, proposed in FigEffa)] initiates the ICD process. 

The starting point of the ICD, in which the two elec¬ 
trons are in the conduction band is presented in KigC |[ l>)| 
One can see that the electron located in the left well is 
in an the excited state while the other electron is in the 
bound state located in the right well. As a result of the 
correlation between the two electrons, i.e. repulsive inter¬ 
action, the electron in the left well decays to the ground 
state of the system transferring its excess energy to the 
electron on the right well which is ionized as a result 
of the electrons correlation (see Fig®. By applying 
low bias to the QWs nano-structure we should sense the 
change in the dark current due to the ICD process. 


ICD process. Second, the electronic correlation has to be 
large enough to make the ICD effective. Since the corre¬ 
lation depends on the distance between the electrons we 
need the wells to be close enough for the correlation to be 
effective. Therefore, on one hand the wells need to be far 
apart to prevent from tunneling to prevail, on the other 
hand the wells should be close enough to allow effective 
correlation between the electrons. The third requirement 
is determined by the conservation of energy in the pro¬ 
cess. Thus we need the relaxation energy of the electron 
in the left well to its ground state to be larger than the 
ionization energy of the electron in the right well. This 
is shown schematically in Figj2]i.e. A El > A Er. 


III. METHODS 



(a) Schematic presentation of the valence and conduction bands of a 
double QW nano-structure composed of the materials shown in 
FigQ] To initiate the ICD process we need to excite the two 
electrons from the valence band to the conduction band. Each 
electron is excited with a different laser frequency uj j, ui 2 



(b)Here we focus on conduction band of the double QW 
nano-structure. Due to the excitation, one electron is in an excited 
state located on the left well and the second electron is in the 
bound state located on the right well. 



(c)Due to the correlation between the electrons, the ICD process 
occurs. The electron in the left well transferring its extra energy to 
the electron in the right well which is ionized. 


FIG. 2 . The I CD proce ss i n th e two coupled quantum wells 
nano-structure (a) (b) and (c) 


The one dimensional effective Hamiltonian we use con¬ 
sists of two electrons in the double quantum well nano¬ 
structure shown schematically in Fig|2] The interaction 
between the two electrons is a soft coulombic repulsion. 
We take into consideration that the electrons in the ex¬ 
periment do not change their momentum in y and z di¬ 
rection. In our calculation we include only the dimen¬ 
sion parallel to the layers growth direction x (see FigQ]), 
in this direction a double QW nano-structure is formed. 
The Hamiltonian of the system is given by : 

H(x i,x 2 ) = H 0 (xi,x 2 ) + V int (x i,x 2 ) 

= h(xi) + h(x 2 ) + V int (xi,x 2 ) (1) 

where the electrons’ positions are represented by xi,x 2 - 
The Hamiltonian in EqQ] contains a sum of two non¬ 
interacting one electron hamiltonians h(xi),h(x 2 ) which 
are coupled by the interaction between the electrons Vi n t ■ 
The one electron hamiltonian is given by : 


h(xi) = P{Xi)- 


eff( X i) 


P{Xi) + V(X,) 


( 2 ) 


The first term in EqO is the kinetic energy term where 
P[xi) is the momentum operator. The kinetic energy op¬ 
erator takes into account the effective mass of the elec¬ 
tron in the different semiconductor layers of the QW de¬ 
picted in FigQ] i.e. this effective mass presented in Eq(2] 
as m e ff{xi) and it is discontinuous in x. The second term 
in EqQ] is the potential energy in the conduction band of 
the double QW nano-structure represented schematically 
in FigQjpy c) This is an effective potential that takes 
into account only one conduction band, it will be given 
explicitly in the next section. The interaction between 
the electrons Vi n t{x\,x 2 ) in EqQ) is the soft coulombic 
repulsion between the electrons and is given by: 


There are few requirement for the ICD process. First, 
the bound state wave-functions of the electrons in the two 
wells should not overlap each other. If the bound states 
overlap then tunneling of an electron from one well to 
the other can take place. This process can overcome the 


Vint (^l 5 x 2) 


_A_ 

y/ (xi — x 2 ) 2 + aexp[—[3(xi — x 2 ) 2 ] 


( 3 ) 
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The interaction strength A in Eq{3] is set from the rela¬ 
tive permittivity of the layers. For the semiconductors 
used in our system (see Fig[T|) the relative permittivity is 
£ r = lOeo, where £o is the vacuum permittivity [dpi. The 
last term in the square root in Eq[3]is introduced to avoid 
the singularity of the potential at x\ = X 2 but to keep 
the coulombic force when the electrons are far apart. The 
ICD process is a decay process which has a finite lifetime. 
We want to calculate the decay rate of this process and 
insure that the lifetime is at list on the same timescale as 
the shortest decay process in the system. We are calcu¬ 
lating the decay rate of the ICD process using the Fermi 
golden rule formula [2^|[27j|. such that it is given by: 

0-7r ~ 

r = = K \('Zf(x 1 ,x 2 )\H\*. l (x 1 ,x 2 ))\ 2 p(E c ) (4) 


The lifetime of the ICD process is given by r = -jt. The 
functions 'F/, Hq are the eigen-state wave functions of the 
unperturbed hamiltonian Hq defined in Eq[T| 4q repre¬ 
sents the initial step of the ICD where the two electrons 
are in the bound states of the corresponding wells (see 
Figl ^b)[ ). T/ depicts the final step of the ICD where one 
electron is in the ground state and the second electron is 
ionized to the continuum (see Fig®. These functions 
are either symmetric (singlet)or anti-synrmetric (triplet) 
with respect to the exchange of the two electrons. These 
are formed from the eigen-state wave functions of the 
one-electron hamiltonian proposed in EqJ5] and are given 
by: 


'S’ i {xi,x 2 ) = ± i>b2{x 2 )i>b 1 (xi)\(5) 

^f(xi,x 2 ) = ^[^ l {x 1 )^ c {x 2 )±^ 1 {X 2 )^ c {x 1 )] 


V’fci i are ^he two bound state wave-functions localized 
in the left well of the double QW shown in Fig J5] with en¬ 
ergies of E^ 1: E^ 2 respectively. V’&i is the ground state 
wave function of the QW, and ipffo is the excited state 
wave function, is the bound state wave-function lo¬ 
calized in the right well with energy of Eijj c is the 
continuum state wave function of the electron which is 
ionized from the right well (see Fig®. The energy of 
this continuum state E c is determined by the conserva¬ 
tion of energy in the ICD process, and is given by: 


77>L TT'L _ Tp TT'R 

&b2 ~ ^bl ~ ~ &bl 


( 6 ) 


In the calculation of the decay rate T we assume that 
only one electron is ionized to the continuum and its mo¬ 
mentum in the y, z direction do not change. Accordingly, 
we consider a one electron and one dimension density of 
continuum states p{E). It is calculated at the energy of 
the ionized ICD electron E c given from the conservation 
of energy in the process depicted in Eq[6l In the next 
section we are going to show the results of the calcula¬ 
tions to the lifetime of the ICD process using the Fermi 
golden rule given in Eq(U 


IV. RESULTS AND DISCUSSION 


The lifetime of the ICD process can be manipulated 
by controlling the physical dimensions of the double QW 
potential shown schematically in Fig(5] The dimensions 
of the different semiconductors layers in the QWs nano¬ 
structure were optimized to increase the yield of the ICD 
process, i.e. to reduce its competition with other decay 
processes existing in the system. In our calculations we 
use the parameters of the conduction bands of the QW 
nano-structure suggested in SecJTTl The only parameter 
varied in our calculation is the distance between the two 
coupled QWs. 

To evaluate the lifetime of the ICD process at a spe¬ 
cific distance between the wells we first calculated the 
eigenvalues and eigen-states of the one electron hamilto¬ 
nian proposed in Eqj2j in our calculation we used only 
one conduction band i.e. single band effective mass ap¬ 
proximation [22]. This effective one dimension potential 
in this hamiltonian is a piecewise potential of the double 
QW nano-structure (see Fig|2]). The potential well depth, 
and the effective mass of the electrons in the different lay¬ 
ers is based on the QW structure proposed in FigfTJ We 
calculated this hamiltonian’s eigenvalues by demanding 
the continuity of the function and flux at the discontinu¬ 
ity points of the potential. The bound states wave func¬ 
tions were required to be square-integrable functions, and 
by using the transfer matrix method we calculated also 
the scattering state wave functions of this hamiltonian 
[361 [37| . To ensure that the effective one dimensional po¬ 
tential is accurate enough, we compared the bound state 
energies and wave functions with a calculation done us¬ 
ing the k ■ p method J28|-j32}]. In this method we take 
into account eight bands from the conduction and valence 
bands, and calculate the electronic bound state energies 
and wave functions. 

The three bound state wave functions and their ener¬ 
gies in the double QW nano-structure are shown in Fig(5] 
for a specific distance between the wells of d = 163A We 
present in Fig[3l the bound state wave functions and en¬ 
ergies from the the two methods. First is the calculation 
from the effective one dimension potential (upper panel), 
this is the bound state wave functions we use to calculate 
the ICD decay rate. Second, to compare our single band 
effective mass approximation, we calculated the bound 
state energies and wave functions using the k • p method 
(lower panel). In both calculations presented in Fig[3]the 
depth of the wells is Vq = —500 mev, while the widths of 
the left and right wells are Wl = 53.96A, Wr = 20.05A 
respectively. 

We obtained in both calculations that the left well sup¬ 
ports two bound states while the right well supports only 
one bound state (see Fig(3|). The ground state wave func¬ 
tion ip^\ and the bound state wave function located in the 
right well in both calculations are almost identical. 
The excited bound state wave function located in the left 
well b’fc 2 Fas a mild difference between the two calcula¬ 
tions. The bound state energies are presented in TableU 
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FIG. 3. (color on line) The one-dimensional effective potential 
presenting the conduction band of the coupled QWs nano¬ 
structure. The depth of the well is Vo = — 500meV, the 
width of the left and right well respectively is Wl = 53.96 
A, Wr = 20.05 A, the distance between the wells i.e. the 
barrier width is d = 163 A (solid red-line).The upper panel 

- shows the calculation of the bound state energies and wave 
functions using the constant piecewise potential i.e. using the 
single band effective mass potential. The lower panel - shows 
the bound states energies and wave functions using the k ■ p 
approximation, which takes into account eight bands from 
the conduction and valence bands. In both calculations the 
effective mass discontinuous in x. The first bound state wave 
function located on the left well ip bl (solid blue).The second 
bound state wave function located on the left well (solid 

- green line). The bound state wave function located on the 
right well i/’bi (solid-black line) 


using the two approaches, E£ b presents the energies cal¬ 
culated using the one dimension effective potential, and 
E^' p (mev) presents the energies calculated using k ■ p 
method. In the last column of Table]]] we present the 
percentage of the conduction band contribution to the 
bound state energies and wave functions x c b us i n g the 
k ■ p method. By comparing the results from both ap¬ 
proaches we obtained that the ground state energies E ^ 
are identical. The energies of the bound states located in 
the right well are very close. In these two bound state 
energies we can see that \cb P — 90 is very high. There 
is a bigger difference in the energy of the excited bound 
state located in the left well Efr 2 , this is due to the fact 
that the percentage of the conduction band contribution 
to that bound state is only = go. 

From this comparison we can see that the bound state 
energies and wave functions are very similar, there is a 
slight shift in the bound state energy containing more 
valence band contribution, this changes will not change 
the result of the ICD decay rate significantly. Further¬ 
more in both methods changing the distance between the 
wells d did not change the bound state energies. This 
comparison between the bound states obtained from the 


Ipbs 

E£ b (mev) 

Efr' p (mev) 

k-p 

Xcb 

Ipbl 

-388 

-388 

0.91 

V’fi 

-210 

-222 

0.88 

V4 

-69 

-128 

0.8 


TABLE I. The three bound state energies using the single 
band effective mass approximation i.e. taking into account 
only one band from the conduction band of the double QW 
(E£ b ). The bound states energies using the k.p approximation 
i.e. taking into account eight bands from the conduction and 
valence bands of the double QW (E bp ). The percentage of 
the conduction band contribution to the bound states using 
the k.p approximation (Xci, P ) 


k ■ p method and the one-dimensional effective piecewise 
potential has shown that we can use the single band ef¬ 
fective mass approximation in the calculation of the ICD 
decay rate. 

In all the calculations of the ICD lifetime all the poten¬ 
tial parameters remained constant and only the distance 
between the wells was varied. We calculated the energy 
of the ionized electron, i.e. the energy of the continuum 
state ip c in Eq[5j The energy of this continuum state 
is E c = 109.3meV which is derived from the conserva¬ 
tion of energy requirement in the ICD process given in 
Eq]6] Although the bound-state energies and wave func¬ 
tions are hardly changed with the distance between the 
wells, the spacial structure of the continuum state wave 
function changes with the distance. In FigQ] we show 
the two continuum state wave functions at two different 
inter-well distances. As one can see the continuum state 
amplitudes and shapes are changing with the distance 
between the wells. These changes affect the ICD life¬ 
time, due to the fact that the overlap between the bound 
state wave functions vanishes and therefore only the con¬ 
tinuum state wave functions couples between the bound 
state wave functions. 


We calculated the decay rate of the ICD process using 
the fermi golden rule formula in EqU while changing the 
distance between the wells. Here we separated between 
the singlet and triplet eigen functions of the unperturbed 
hamiltonian in Eq]5] In Fig 5(a) we show the lifetime 
of the ICD as a function of the distance between the 
wells using the triplet functions, we got similar results 
also for the singlet functions. We expect the life time 
to grow as the distance between the wells increases, due 
to the fact that as the distance increases the correlation 
between the elect rons decreases and so does the decay 
rate. In Fig |5(a)| we see that the overall trend follows 
this expectation except around the point of d = do = 
136. 8A. At this point we see a surprising sharp drop in 
the lifetime. One can see that although the distance at 
this point is quite large, we get very short life time of 
several picosecond. This is an order of magnitude shorter 
then what is expected. 

To show that this is indeed a result of ICD we also cal¬ 
culated the overlap between the bound state wave func- 
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FIG. 4. (color on line) The continuum state wave functions 
of the ionized electron ip c due to the ICD process. The differ¬ 
ent panels referring to different distances between the wells in 
the one electron potential given in Fig[3] upper panel - the 
distance between the wells is d = 163A, lower panel- the dis¬ 
tance between the wells is d = 136.8A. The continuum state 
has energy of E c = 109.3meV in both panels, this energy is 
set from the conservation of energy in the ICD process see 
Eq[6] One of the continuum state ip c describes the electron 
arriving from oo while the other continuum state describes 
the electron that arrives from — oo (solid black and dashed 
magenta respectively). The edges of the left and right wells 
in the double QW nano-structure are shown in solid blue and 
solid red lines respectively. 


tions of the right well with those on the left well by treat¬ 
ing th e well s as separate systems. The results are shown 
in Fig |5(b)| As one might intuitively expect, we see that 
as the distance between the wells increases the overlap 
between the bound state wave functions in the two dif¬ 
ferent wells decreases. This overlap is a measure of the 
tunneling in the system, which means that tunneling is 
effective at small distances. When tunneling occurs in 
our system the ICD is not the dominant decay process. 
Around the special point of do the overlap between the 
bound state wave functions of the two wells is very small 
so we know that the tunneling is not effective. This is in 
contrast to smaller distances where the lifetime is short 
but the overlap is large enough to make tunneling the 
dominant process. 

In order to explain this interesting result, we need to 
examine the expression for the ICD decay rate given in 
EqQ] One possible reason for the increased decay rate 
could lie in the shape of the continuum wave function 
in the region between the wells. We can look at the 
continuum state amplitude of two different distances: at 
do and at d\ = 163A (see FigHJ), and observe only a small 
difference in the amplitudes of these two wave functions. 
This cannot be the reason for such a significant change 



(a)The lifetime of the ICD process r as a function of the distances 
between the wells d. 



20 40 60 80 100 120 140 160 180 

d (A) 


(b)The overlaps between bound state wave functions of the two 
wells (treating each well separately) as a function of the distance 
between the wells. The overlap of the bound state in the right well 
and the ground state in the left well (solid red). The overlap of the 
bound state in the right well and the first excited state in the left 
well (dashed blue) 

FIG. 5. The ICD life time is shown in Fig |(a)[ and the overlap 
between the bound state wave functions in the different wells 
is shown in Fig|(b)| 


in the lifetime. 

Another factor that may affect the decay rate is the 
density of states in the continuum at the energy of the 
escaping electron p{E c ). The density of states in the free 
one-electron picture follows p/ ree oc If this was 

the case in our system the trend in the lifetime should 
not change because the energy of the escaping electron 
E c = 109.3meV remains constant with the distance. Due 
to the structure of the potential, the density of states 
has peaks which are correlated with this structure. We 
calculated the density of continuum states for the one- 
electron hamiltonian in Eqj2] (used also for the lifetime 
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calculations) at different distances by choosing vanishing 
boundary condition in a large box ijj(x = ±L/2) = 0 
[39l [40 ]. We evaluate numerically the one dimensional 
density of states in the energy of the continuum states i.e. 
energies above the threshold, by calculated the following 
expression in a large box: 


p{E) 


-(—r 1 

L V An’ 


(7) 


We normalized this density of states by dividing it to the 
length of the box L. 

Fig | 6 (a)| shows the density of the continuum states of 
the one-electron hamiltonian h(xi) as a function of the 
continuum energy in two different structures, i.e. two 
distances between the wells. The solid blue line is at the 
distance do , and the red dashed line is at a distance of 
di = 163A. One can see that although the background 
remains the same in both plots, the peaks appear at dif¬ 
ferent energies. Furthermore at do there is a peak in the 
density of states at the energy of the ionized IC'D electron 
E c . While in the other distances such as in d\ the energy 
of the ionized ICD electron which is also E c there is no 
peak in the density of states. This means that the decay 
rate will be an order of magnitude larger than expected 
in the spacial distance of do- This explains the sharp 
drop in the ICD lifetime in Fig |5(a)| around the point do- 

One can calculate the resonance states of the one- 
electron hamiltonian h{xi) given in Eqj2]by imposing out¬ 
going boundary conditions on the schrodinger equation 
[381 ]. The resonance picture of the one electron hamil¬ 
tonian can help in understanding the density of states 
picture. The resonance energies on the complex energy 
plane for inter-well distance of do are shown in Fig | 6 (b)| 
It is evident that the position of the resonance ener gies 
matches the peaks in the density of states (see FigEfail 
in solid blue). 

This can be explained by the fact that the resonance 
states are connected directly with the poles of the S- 
matrix in the complex plane |38|. When the poles of the 
S-matrix are isolated from each other, and close enough 
to the real axis, one can associate the peaks in the cross 
section with real part of the poles of the S-matrix on the 
complex energy plane [4lj . The cross section depends on 
the density of states via the S-matrix or Green operator 
[i^ ]. such that the peaks in the density of states should 
appear at the energies of the poles, i.e. the positions of 
the resonance states. 

From this discussion we conclude that the efficiency of 
the ICD process is enhanced at the distance of do- This 
makes the ICD life-time on the same time scale of the 
dominenet competing decay process in the QWs nano¬ 
structure. The reason for this significant enhancement is 
that the ionized ICD electron is temporarily trapped in a 
shape type resonance state. This trapping enables us to 
get an efficient ICD ionization that competes with other 
decay processes even at very long distances. Therefore 
the calculations of the shape type resonances as function 
of the distance between the two QWs provides us a pow- 



(a)The density of continuum states for the one electron hamiltonian 
h(xi) of the system for two different inter-well distances, calculated 
numerically by solving the problem in a large box L = ! 5870/1. The 
density in a distance of d = do = 136.S/1, which is the spacial 
distance where there is a sharp drop in the lifetime of the ICD (see 
Fig|5( a)J (solid blue). The density for a distance of d\ = 163A 
(dashed red). We can see that that the background of this plot is 
the density of free particle in one dimension. The peaks in the 
density are correlated with the structure of the potential, therefore 
at different distances we see different picture of the density of 
continuum states and its peaks. The continuum state energy of the 
ICD electron escaping the system is marked on the plot as E c and 
it do not change with the distance. 



Re (E) 


(b)Resonance solutions of the one electron Hamiltonian in Eq[2on 
the complex plane with a distance of d = do = 136. 8v4 between the 
wells. One can see the correlation between the resonance position 
and the peaks in the density of states 


FIG. 6. (a) the density of continuum states for the one elec¬ 
tron hamiltonian presented in Eq[2]m two different inter-well 
distances, (b) The resonance solutions to the one electron 
hamiltonian on the complex plain. 


erful computational tool for designing an experiment to 
observe the ICD phenomenon. 

























V. CONCLUDING REMARKS 

In this work we propose for the first time an ICD pro¬ 
cess in a realistic structure of two coupled QWs. The 
system we studied is based on real physical parameters 
of QWs in semiconductors materials. The one electron 
effective potential is set from the band-structure of the 
semiconductor materials in the QW nano-structure. We 
also take into account the discontinuity in the effective 
mass of the electrons in the different QW layers, and 
the permittivity of the layers. The ICD lifetime depends 
very strongly on the distance between the wells. The 
overall trend is that by increasing the distance between 
the wells the lifetime is increasing. Our results shows 
that the shortest ICD lifetime in our system is of several 
picoseconds. This means that the ICD is on the same 
time scale of inter-subband relaxation with LO phonons 
which dominant inter-subband relaxation process in QW. 
By designing a sample which matches the ICD conditions 
and with enough double QW periods, this phenomenon 
should be observable experimentally. 


The main result of this paper shows that the param¬ 
eters of the potential can be manipulated such that the 
ICD process is enhanced. The ionized electron is tem¬ 
porarily trapped in a shape type resonance state, this 
resonance state introduces a peak in the density of con¬ 
tinuum states resulting in a very short lifetime of the 
ICD process on the time scale of picosecond, even at a 
very long distance between the wells. Based on our re¬ 
sult and understanding of the ICD process, we can design 
an experiment which will show this phenomenon for the 
first time in nano-structures. This can lead to designing 
a photo-detector which is very sensitive to wave-length 
even at very low intensities . 
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